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^ ■ Abstract 

A method for calculating the retarded Green's function for the gravitational wave equation 
in Friedmann-Roberson- Walker spacetimes, within the formalism of linearized Einstein 
gravity is developed. Hadamard's general solution to Cauchy's problem for second-order, 
linear partial differential equations is applied to the FRW gravitational wave equation. 
The retarded Green's function may be calculated for any FRW spacetime, with curved or 
flat spatial sections, for which the functional form of the Ricci scalar curvature R is known. 
The retarded Green's function for gravitational waves propagating through a cosmological 
fluid composed of both radiation and dust is calculated analytically for the first time. It 
is also shown that for all FRW spacetimes in which the Ricci scalar curvatures does not 
vanish, R ^ 0, the Green's function violates Huygens' principle; the Green's function has 
support inside the light-cone due to the scatter of gravitational waves off the background 
curvature. 







I. Introduction 

The Green's function is an important tool for evaluating the solution to a partial 
differential equation for a given set of boundary conditions. The Green's function repre- 
sents the solution to Cauchy's problem: given a set of properly formulated initial data, the 
Green's function may be used to uniquely determine the final, evolved data. In physics, 
one may determine the evolution of a field according to its equations of motion, given some 
initial data. The techniques used in this paper have been widely applied to the calculation 
of Green's functions, or propagators, for quantum fields in curved spacetime. Here, those 
techniques will be applied to classical fields in cosmology. In this paper I will determine the 
evolution of a tensor field, representing gravitational waves, given a source stress-energy 
tensor. The goal of this paper, then, is to calculate the Green's function for gravitational 
waves evolving in an arbitrary FRW spacetime. 

The Green's functions may be useful for the study of gravitational radiation in FRW 
spacetimes, especially as it relates to observational cosmology (for an example, see [1]). For 
situations in which the stress-energy tensor of a source of gravitational waves is known, the 
Green's function may be used to calculate the properties of the tensor metric perturbations, 
and the spectrum of gravitational radiation. Examples of such sources are cosmic strings 
or scalar fields present in the early universe. The tensor perturbations generated by such 
sources may be manifest as anisotropies in the cosmic microwave background, or as a 
spectrum of stochastic gravitational radiation. Then, the Green's function may be used 
for the calculation of physically observable quantities. 

The organization of the paper is as follows. In section II, the formalism for calculating 
the Green's function, following Hadamard's general solution to Cauchy's problem for the 
wave equation, will be presented. A similar technique was developed by Waylen in [2,3]. In 
section III, this technique will be applied to find the Green's functions for arbitrary FRW 
spacetimes. The main result of this paper, a general expression for the Green's function for 
any FRW spacetime in which R is known, will be presented. Several specific cases will be 
evaluated explicitly, including the case in which the cosmological fluid undergoes a smooth 
transition from radiation- to dust-dominated expansion. In section IV the properties and 
applications of these Green's functions will be discussed. The focus will be on the cosmo- 
logical significance of those terms which violate Huygens' principle. A final summary of 
this work will be presented in section V. 

II. Formalism for calculating Green's functions 

The technique for calculating Green's functions follows from Hadamard's general so- 
lution to Cauchy's problem for a second-order, linear partial differential equation [4]. The 
most general form for the solution to the tensor wave equation 

h ab . c c + Ah ab = (III) 

may be written as an expansion in powers of the non-local, biscalar of geodetic interval 
a(x,Xi) = ^s(x,Xi) 2 , where s(x,Xi) is the geodetic interval between a fixed 4-vector x^ 
and a free 4-vector x in the spacetime. (Note that A is some arbitrary function, not the 
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cosmological constant.) The Green's function solution is a bitensor, an object with tensor 
indices on the spatial hypersurfaces Xi and x: 



G(x, Xi ) ab Cidi = ^ [u ab c ^ i + v ab c ^ log \a\ + w ab c ^} . (II. 2) 

Here, u, v, and w are bitensors which are free of singularities. Extending II.2 into the 
complex-a plane, and retaining only the imaginary part, we may obtain the retarded 
solution [5,6] 

G ret (x, Xl ) ab c ^ = ±-[u ab c ^5{<j)-v ab c ^e(-<j)\e(\). (11.3) 

The parameter A runs between the spacetime points Xi and x, such as the time, and is 
positive when the spacelike hypersurface containing x lies to the future of that containing 
Xi. In that this equation represents the solution to the wave equation which is homogeneous 
everywhere except at a single point, we may obtain the solution to the inhomogeneous wave 
equation 

h ab -/ + Ah ab = f ab (II A) 

by summing the contribution of many points: 

h(x) ab = [ ^(^)d 4 x'G ret (x,x') ab Cidi f(x') Cidi . (II.5) 

Thus, the Green's function may be used to construct solutions to the wave equation. 
Equation i7.3, then, gives the retarded Green's function which we seek. 

A prescription for calculating the non-singular, bitensor quantities u and v may be 
obtained by generalizing the work of DeWitt and Brehme [5,6] and McLenaghan [7] on the 
scalar and vector wave equations. Such a calculation has been carried out by Waylen [2,3]. 
One may apply Hadamard's general solution II. 2 to the homogenous wave equation, and 
equate terms with matching powers of a. One obtains 



(a. c c -4)u ab c ^+2u ab c ^, c a^ = 
2a. c v ab c ^' c + (a, c c - 2)v ab Cid > + u ab c ^ ./ + Au ab c ^ = 0. 



(JJ.6) 



Manipulating these equations, we find that [7] 
lnu = - l -( d 4(v./-4) 



A 

Va^ = - ^U ab ^ J d\(u ef ^./ + AlX^XO 6 ' 9ihi 



(II.7) 



In general, the bitensors u and v will be proportional to the bitensor of parallel geodetic 
transport (see the discussion of such bitensors in [5]). The functions u and v may be 
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evaluated once the biscalar of geodetic interval a is known. Therefore, for a given metric, 
a, u, and v may be determined, and the Green's function for the tensor wave equation 
may be evaluated. This general technique may be applied to the wave equation in any 
spacetime. 



III. Green's functions in FRW spacetimes 

The gravitational wave equation for which the Green's function will be determined, 
describing the propagation of tensor metric perturbations on a background spacetime, may 
be derived in the context of linearized, Einstein gravity. (For some details on linearized, 
Einstein gravity, see [8,9]). I will consider tensor metric perturbations, g ab — > g ab + h ab , 
where the line element for the background, FRW metric is 

ds 2 = a( V ) 2 (-dri 2 + y~k^ + r2d ° 2 + r ' sin2 ed ^- 

Due to the symmetries of the FRW spacetime, I may apply the transverse, traceless, 
synchronous (TTS) gauge conditions [10,11] 

h ab ' b = h ab g ab = h ab t b = (III.2) 

to isolate the physical degrees of freedom in the gravitational radiation. Here, t a is a unit 4- 
vector pointing in the time direction ?y, orthogonal to spacelike hyper surfaces. Ultimately, 
the gravitational wave equation is found to be 

16nGT:r ce = h ab ;c c + 2R C a d b h cd . (III.?,) 

Here, R c a d b is the Riemann curvature tensor in the background FRW spacetime, and G is 
the gravitational contant. (We use units such that the speed of light is unity: c = 1.) The 
stress-energy tensor T^ urce represents the external source of the gravitational waves. 

We may now apply the techniques developed in the previous section in order to eval- 
uate the retarded Green's function. The bitensor functions u and v are proportional to 
the transverse, traceless, synchronous bitensor E(x, Xi) ab Cid \ which projects only the TTS 
portion of the source stress-energy tensor from the spacelike hypersurface at x, to x. The 
TTS bitensor T, ab Cidi may be constructed through the use of the spatial biscalar of geodetic 
interval, fi(x,Xi) = \x — Xi\, following the work of Allen and Jacobson [12,13]. Defining 

V a = V a At(x , X i) Pab=9a b + t a t b < . 

{Ill A) 

V ai = V ai v(x,Xi) P ahi = g b .P ab 
the unique TTS bitensor with indices a, b on the tangent space at x, and q, di at Xi is 

Z ab ° idi =P a C *P b di + P/'Pb ' ~ P ab P Cidi + V a V b P Cidi + P ab V C *V d > 

+ P a Cl V b V di + P a di V b V c > + P b Cl V a V d * + P b dl V a V Ci + V a V b V c *V di . 
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Here, V a E ab Cidi = g ab E ab Cldi = t a Y> ah Cidi = 0, such that it is transverse, traceless, and 
synchronous. Thus, combining equations II. 3, 7, and III. 3 the entire prescription for 
calculating the Green's function for tensor metric perturbations is specified. 

The calculation of the Green's function for FRW spacetimes may now be carried 
out. Having specified the background metric, the biscalar of geodetic interval is a = 
\(\x — Xi\ 2 — (rj — i]i) 2 ). Evaluating equations 77.7, 

cidi Cidi 
almjalri) 



12(7/ - Tji) J n 

In this equation C(x) = xj sin(x) and 1Z is the curvature scale of the spatial sections of 
the spacetime such that S R = 6TZ~ 2 . In the case K = 0, 1Z — > oo so that C = 1. Applying 
the above results to equations II.3 and JiT.3, one may write the Green's function as 

t( , c% d % C{\x-x'\/n)\5{'q-r] i -\x-x i 
G {x,Xi) ab =- 



47ra(i]i)a(r)) 



+^ ( f a 2 ( v )Rd v ) e{v m ls M l eh - thWx, Xi ) ab c ^ . {IIL7) 

Note that the Green's function is symmetric under interchange would be 

expected. Also, the Green's function possesses the correct dimensionality, (length) -2 , 
where the coordinates (77, x) are dimensionless and the expansion scale factor a carries 
units of length. This equation composes the main result of this paper. 

The Green's function may be evaluated for a number of specific cases. The results are 
displayed in the accompanying table. (See table.) In this table, a = —1, 1,2 for deSitter 
spacetime, radiation- and dust-dominated expansion, respectively. Here, the equation of 
state of the cosmological fluid is p = ^^P- Minkowski spacetime is given by K = 
and a(rj) = 1 (for which case the preceding equation of state is invalid). The first three 
entries in the table give the expansion scale factor and the Green's function for an ideal 
cosmological fluid, in a spacetime with flat and curved spatial sections. The results from 
the first entry, for flat spatial sections, have been previously calculated for the specific 
cases a = —1,1,2 [1,14]. The Green's function for arbitrary a and for any K is given 
here for the first time. The last three entries in the table give the expansion scale factor 
and Green's function for a mixed cosmological fluid, in a spacetime with flat and curved 
spatial sections. In the case of the mixed ideal fluid, the energy density is p = p ra d + Pmat 
and the pressure is p = \p r ad, where r is a function of p cr it(vo) an d p r ad(vo)i an d rjo is 
the present-day conformal time. Then, for ij<t, the fluid is radiation-dominated. For 
Tj 3> t, the fluid is dust-dominated. The Green's function for this mixed case is presented 
here for the first time. This case represents the most realistic cosmological scenario. 
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IV. Properties and applications of Green's functions 

The Green's functions describing the evolution of gravitational radiation in FRW 
spacetimes, derived in the preceding section, display interesting properties and have useful 
applications, which may now be discussed. Specifically, one may examine the physical 
consequences of the terms in the Green's functions which contribute to the violation of 
Huygens' principle. Applications of the analytic Green's functions for open and closed 
spacetimes, for a combined radiation and matter fluid will be discussed. 

Huygens' principle states that the effect of a luminous disturbance about some point 
at time t will be localized at a later time t' in a very thin spherical shell at a radius c{t' — t) 
[4]. So, a wave packet travels spherically outward on the light cone from some source. 
From examining equation II. 3, one finds that for v 7^ 0, the Green's function has support 
inside the light cone; the "wave packet" disperses as it travels and is not localized on a thin 
spherical shell. In an FRW spacetime, equation III. 7 indicates that when R 7^ Huygens' 
principle is violated. 

The violation of Huygens' principle is due to the coupling of the curvature of spacetime 
to the tensor metric perturbations. The energy in the background curvature feeds into the 
gravitational radiation. This coupling of spacetime curvature to fields has been well studied 
on the microscopic level in the context of particle creation [15,16]. (See references [17,18] 
for a discussion of quantum mechanical propagators with respect to Huygens' principle.) 
A classical analogue of particle creation is observed here on the macroscopic level [19,20]. 
Examining the table, the function V increases with the expansion a(rf); V represents the 
amplification of the classical gravitational field due to the expansion. One may say that 
energy is transferred from the expansion to the gravitational field. Ultimately, it may 
be simply stated that Huygens' principle for gravitational waves propagating in an FRW 
spacetime is satisfied only when R = 0, or the spacetime is filled with a conformally- 
invariant radiation fluid. 

The Green's function for the case of the mixed radiation-plus-dust fluid, for K = 
±1,0, as presented in the table, is new. These analytic expressions should be useful for 
examining the evolution of gravitational waves through the transition from radiation- to 
dust-dominated expansion. 

The primary application of the Green's functions derived in this paper is for the 
situation in which the source stress-energy tensor is known. The source ought to be "stiff" , 
in that the source evolves freely of the background spacetime and the perturbations it 
produces. Some examples are cosmic strings, global topological defects, and scalar fields 
[1,14,21,22]. Of interest are the anisotropies produced in the microwave background, and 
the spectrum of gravitational radiation emitted. 

The anisotropies in the microwave background may be calculated by summing the con- 
tributions to the temperature fluctuations caused by the perturbations in the gravitational 
field along the path length of a photon traveling from the surface of last scattering. 

AT 1 f 1 

= -\J dXea ^ eb ^^{J(^ J ^id A x i G™\x,x i ) a S idi l^GT s ° urc *{x i ) Cid )) 

(IV.l) 
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For a given source stress-energy tensor, the temperature fluctuations 5T/T causedd by 
cosmic strings, global topological defects, or scalar fields [1,21,14,22] may be calculated. 

The energy density and power in the gravitational radiation emitted by a given source 
may be examined by calculating the stress-energy tensor 

8llGT^ aV =-h cd [hda;bc + hdb;ac ~ h ab - cd — h c d- a b] ~ d; a h d C ;b 

1 1 K ' 2 {IV - 2) 

+ -9ab [h Cd h cd , n n + -h cd ; n h Cd -> n - h Cd h cd {- + \)] • 

2 L 4 a z a 4 

This stress-energy tensor, the terms occuring in the perturbed Einstein's equations which 
are second-order in h, may be evaluated for a given source through the use of the Green's 
function. This expression may have use in studying the power in tensor metric perturba- 
tions produced during inflation (for a recent review, see [23]), or through the collision of 
bubbles formed in a first-order phase transition [24]. In this case, the Green's functions 
derived in the preceding section for the mixed ideal fluid may be especially useful for cal- 
culating the evolution of the gravitational waves through the transition from radiation- to 
dust-dominated expansion. 

Finally, it is interesting to note that while the Green's function for the mixed ideal 
fluid has been obtained, the mode functions for the gravitational waves in such a case 
apparently cannot be written in terms of known functions. In the case of a single, ideal 
fluid background, the mode functions, solutions to equation iT.l, take the form of Bessel 
functions. The case of the mixed fluid, however, appears much more complicated. 



V. Conclusion 

In this paper the retarded Green's functions for the gravitational wave equation in 
Friedmann- Robertson- Walker spacetimes, within the formalism of linearized Einstein grav- 
ity, were calculated. While the form for the Green's function for a generic FRW spacetime 
was presented, several specific cases were considered. These cases, presented for the first 
time, include an ideal cosmo logical fluid with an equation of state p = ^^p, and the 
case of a mixed ideal fluid of radiation and collisionless dust, where p = p rad + p ma t and 
p = ^Prad- The Green's functions for varying spatial curvature, K = 0, ±1 were also con- 
sidered. It was also shown that for all non-conformally invariant FRW spacetimes, in which 
R 7^ 0, the Green's function violates Huygens' principle. This is the classical analogue of 
particle creation in a varying gravitational field, as the gravitational waves scatter off the 
background curvature and gain energy from the cosmological expansion. Finally, it was 
indicated how these Green's functions may be applied to the calculation of the microwave 
anisotropies and spectrum of gravitational radiation produced by "stiff" sources. These 
applications of the Green's functions will be carried out in a future work. 
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Table 



Green's functions for gravitational waves in FRW spacetimes 



G ret {x,Xi)ab 



Cidi _ C(\x-Xj\/1Z) Y S{r)--q l -\x-x l \) 
4-ira(rn)a(r]) L \x— Xi\ 



+ 



y 6(jj-ru-\x-Xi\) 
V-Vi 



]0(r] - r]i)T,(x,Xi) ab c 



K 



Cdx-Xil/ft) 



a(r/) 



cosmology 



a(a-l)(v-Vi) 



ideal fluid 



^/sinh(^) 



/ u sinh(T7/a) i a 
"WvLsinh^i/a)] 



Kn, — 1} ( sinh[(T7-T ?1 )/q] \ 
2^" - L Asinh(? ? /a)sinh(?7i/a)y' 



ideal fluid 



' — Xi \ I ■ I \X — Xi \ 



( Nr sin(77/a) i 



±(rv — 1} ( sin [( ? ?- ? ?»)/ Q ] \ 
2^ - L Asin(? 7 /a)sin(?7 l /Qf)^ 



ideal fluid 



1 l oe - v(Vi+r) 
t ° Vi(v+ T ) 



mixed fluid 



n 



7sinh(l^i) 



/ \ r rsinh?7+sinh 2 (?7/2) i 
" V'V Lrsinh ?7 l +sinh 2 (?7 l /2)J 



1 l n r-r 2rcoth (^/ 2 ) + 1 l 
4r 1U &L 2rcoth(?7/2) + l J 



mixed fluid 



' — Xi \ I • I X — Xi \ 

■^/ sm (V) 



U rsin?7+sin 2 (Tj/2) i 
"VW Lr sin ^+81^(^/2) J 



J_ 1n ^r 2rcotfa/2) + l l 
4r 1U &L 2rcot(r;/2) + l J 



mixed fluid 
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